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ABSTRACT 


An  often  ignored  aspect  of  electromagnetic  radiation  iron 
antennas  is  the  characterisation  of  their  near-fieids.  A 
computer  program,  Numerical  Electromagnetics  Code  (NEC);  is 
validated  for  accurate  near-field  contour  at  ions  and  applied  to 
a  model  of  a  broadcast  mcnop-cle.  E-  and  H-fields  are  plotted 


as  a  sanction  or  position  along  the  antenna  for  various 


distances  from  the  surface. 


’he  fields  are  also  clotted  as 


:unct ion 


idial  distr 
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I.  INTRODUCTION 


A.  PROBLEM  ENVIRONMENT 

Maxwell's  equations  describe,  for  all  space,  the  Electric 
(E)  ana  Magnetic  (H)  Fields  that  are  the  medium  for  trans¬ 
mitting  information  (virtually)  instantaneously  over 
extended  distances.  This  practical  application  of  electro¬ 
magnetic  (EM)  radiation  has  generated  considerable  interest 
in  its  far-field*  characteristics.  The  E  and  H  vectors  can 
be  analytically  determined  when  the  current  (I)  on  the 
radiator  (from  which  the  field  vectors  are  generated)  has 
been  described.  For  the  far-fiald  region  of  monopoles  and 
dipoles,  acceptable  results  are  obtained  by  assuming  a 
sinusoidally  varying  current.  However,  when  the  field 
vectors  are  calculated  for  field  points  at  distances  less 
~han  X  /  2t  (the  near-field  region),  significant  errors  occur 
between  calculated  values  and  empirically  measured  values. 
This  discrecancv  is  not  the  resul*  of  inaccuracies  in 
.Maxwell's  Equations,  but  can  be  wholly  attributed  to  an 
inaccurate  descr ic~ i or.  of  tha  rurrent. 
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In  the  past,  the  inaccurac' as  of  the  near  field  calcula¬ 
tions  were  not  a  problem — there  was  no  practical  application 
to  which  the  energy  could  be  apDlied.  However,  in  recent 
years  there  has  been  an  escalating  interest  in  the  close-in 
problem.  This  interest  has  been  generated  by  an  increasing 
proliferation  of  devices  which  use  Radio  Frequency  (RF)  energy 
in  such  a  manner  that  the  user  is  exposed  to  near-field 
radiation  (e.g.  microwave  ovens,  hand-held  walky-taikies , 
wireless  telephones,  RF  sealers,  etc.). 

The  U.3.  Navy  was  one  of  the  first  organisations  to 
express  an  interest  in  these  fields.  The  interest  was  the 
result  of  a  concern  for  the  hazards  these  fields  might  pcse 
to  fuels,  ordnance,  other  El!  operational  and  test  equipment, 
and  to  personnel.  The  hazards  posed  by  fuels,  ordnance, 
and  other  EM  equipment  result  in  overt  effects;  the  hazards 
posed  to  oerscnnel  are  not  so  obvious. 

With  virtually  unanimous  consensus,  the  literature 
identifies  the  heating  of  boay  tissue  as  the  primary  hazard 
to  personnel.  However,  some  distinctly  ncn-thermal  effects 
have  been  observed.  Rene  of  these  effects  [Ref.  1j  include: 

a)  minor  changes  in  blood  properties 

b)  a  ''buzz"  heard  by  certain  oeople  when  exposed  to 
microwave  radiation* 

‘■•The  "buez"  is  assumed  to  be  a  function  of  the  tulse  repetition 
rate  (RRR)  rather  than  the  carrier  frequenev. 


c)  abnormalities  of  the  chromosome  structure 

d)  movement,  orientation,  and  polarization  of  protein 
molecules 

e)  changes  in  the  transport  rate  at  the  blood-brain 
barrier 

f)  comfort  imbalances,  such  as  epigastric  distress, 
emotional  upset,  and  nausea. 

The  significance  of  these  effects  is  not  yet  understood. 

The  general  public  has  not  been  exposed  to  near-field 
radiation  hazards  of  sufficient  intensify  or  duration  to  be 
cause  for  alarm;  however,  there  are  certain  vocations  in 
which  the  incumbents  have  an  elevated  probability  of 
exposure  to  this  type  of  field.* 

3.  THESIS  STATEMENT,  SCOPE  AMD  LIMITATIONS 

In  order  to  provide  more  accurate  data  for  pecole 
researching  biological  effects  of  near-field  radiation,  this 
study  will  validate  a  computer  program  which  or routes  the 
near  electric,  magnetic,  and  comocsite  (peak)  fields  cf  an 
arbitrary  radiator.  The  program  will  also  be  used  to 
compute  these  near  fields  for  a  model  to  be  described  below. 
The  program,  Numerical  Electromagnetics  Oode  (NEC),  was 
developed  ad  Lawrence  Livermore  Laboratory  under  the 

:!CaiL:rs  at  sea  or  steeole4  acks  who  effect  antenna  remains, 
f” ■> >~>  instance 
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sponsorsmp  or  The  Naval  Ocean  Systems  Center  (NOSC)  and  the 
Air  Force  Weapons  Laboratory.  It  can  be  easily  learned  and 
applied  directly  by  a  researcher  who  has  some  knowledge  of 
electromagnetic  theory. 

The  model  to  be  used  for  the  investigation  will  be  a  75 
meter  high  broadcast  monopcle  radiating  above  a  perfectly 
conducting  ground  plane.  The  ground  plane  is  located  in  the 
X-Y  plane;  the  monopcle  is  co-axial  with  the  F-axis.  The 
antenna  will  be  excited  at  its  base  with  a  voltage  sufficien 
to  "reduce  a  radiated  cower  of  1000  watts.  The  excitation  i 


at  a  frequency  of  1  Megahertz,  which  has  a  wavelength 


or  out 


meters;  hence,  the  monopole  is  a  quarter-wavelength  radiator 
operating  at  resonance.  A  diagram  of  the  model  is  shown  in 


In  simulating  this  model,  NFC  requires  that  the  antenna 
be  broken  into  virtual  shert  straight  segments.  In 
consonance  with  this  requirement,  three  different  configura¬ 
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Field  values  will  also  be  calculated  along  the  radial  axis 
( Y-axis)  in  increments  of  .5  meters  for  heights  (h)  cf  2h, 
ih,  l/2h,  l/4h,  and  2/100ths  h.* 

The  near-field  examined  in  this  study  is  the  cylindrical 
volume  bounded  by  a  cylinder  of  one  wavelength  radius 
centered  on  the  Z-axis,  rising  from  the  ground  plane  to  a 
height  twice  that  of  the  mono cole. 
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II.  DESCRIPTION  OF  SELECTED  ANALYSES 


.  CLASSICAL  ANALYSIS 

The  antenna  utilized  in  this  study  is  a  member  of  the 
lass  identified  as  'thin-wire'  antennas.  The  criteria 
hic.n  adumbrate  the  thin-wire  approximation  are  [Ref.  2]: 

a)  Transverse  currents  are  negligible  relative  to  axia 
currents 

b)  Circumferential  variations  of  the  axial  current  are 
negligible 

c)  The  current  can  be  represented  by  a  filament  currer. 
on  the  wire  axis 

d)  Boundary  conditions  on  the  electric  field  need  be 
enforced  only  in  the  axial  direction. 

By  far,  the  largest  subset  of  thin-wire  antennas  is  a 
lass  known  as  short  antennas.  Chert  antennas  are  define 
s  having  one  dimension  (h)  much  greater  than  the  other 
intensions,  but  h  is  also  much  lass  than  a  wavelength , 
vpicaily  h  <  \ / 3 .  Because  of  the  predominance  of  this 
uc-set,  it  is  reasonable  “hat  it  should  be  us--d  is  a  mod 
or  classical  analysis. 

Based  on  Maxwells  Equations  and  geometric  considerati 
s  shown  in  Figure  2,  the  following  non-cerc  complex  seal 
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where:  r  =  the  intrinsic  impedance  of  the  medium,  I  =  antenna 
current,  X  -  wavelength,  3  =  2v/\  =  phase  constant,  h  = 
monopcle  height,  r  -  distance  of  field  observation  point 
from  monopcle  base.  since  there  is  no  variation  in  f: 


a;  The  distance  to  the  field  observation  point  is  large 
relative  to  the  monopcle  height,  and 

b)  An  harmonically  time  varying  current  is  uniform  over 
the  height  of  the  monopole. 

The  first  of  these  assumptions  is  made  to  focus  attention  to 
that  volume  of  space  where  "he  lines  of  force  have  detached 

“ h*2I?i3  3  1  V3  3  TT  .iTH  DhS  d”  t  find  3T3  3.3  DTCa.Cfli.nS  CH6  5 C TT7!  05 
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a  plane  wave  front.  The  second  assumption  is  to  form  a 
first  approximation  of  the  current  distribution  on  the  wire 
which  is  solvable  in  closed  form.  Extensive  analysis  over 
many  years  has  shown  this  to  be  a  valid  assumption  when  the 
antenna  is  excited  at  a  single  point. 

In  this  study,  we  will  be  interested  in  fields  where  the 
observation  points  are  not  at  distances  considered  large 
relative  to  the  mor.opoie  height.  When  this  is  the  case,  the 
diagram  of  Figure  2  must  be  modified  as  shewn  in  Figure  3*. 
From  the  geometry  of  this  figure,  the  following  non-zero 
complex  scalar  field  components  can  be  derived: 

-  KcosSh)  - - )  (3) 

v* 

*  3 
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0 


( "  > 


* Ecu at ion s  (1)  through  (3)  ere  expressed  in  spherical  coor¬ 
dinates  while  Equations  (3)  through  (')  have  been  transforms 
to  cvlir.drioai  coordinates . 
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Ima of  Moropel 
Reflected  throi/^h 
Ground  Plane 


A  current  or  the  form: 


1(c)  =  I^sinB'n  for  0  <_  z  <_  h 


has  been  inserted  into  Equations  (5)  through  (7). 

Equations  (5)  through  (7)  indicate  that  each  field 
component  is  the  phasor  sum  of  contributions  from  three  coin 
sources.  The  first  source  is  located  at  the  top  of  the  mono 
pole,  a  second  image  source  is  located  in  a  symmetric  pcsiti 
when  the  first  source  is  reflected  through  the  ground  plane. 
The  third  source  is  located  at  the  base  of  the  monopole;  its 
relative  amplitude  is  modified  by  the  (-2cosSh)  factor. 

Since  we  have  a  solution  for  computing  the  field  points, 
it  is  now  reasonable  to  determine  the  range  ever  w.oich  it  is 
valid.  3y  inspection  we  note  that  for  r  =  \/2n  in  Equations 
(1)  through  (;),  the  magnitude  of  each  nen-oero  term  in  the 
oarenthesis  is  equal.  This  implies  that  r  =  \/'2-t  may  be  a 
transition  point.  Further  investigation  indicates  that  for: 

a)  r> \ / 2t :  the  1/r"  terms  in  E„  and  H. ,  and  the  1/r" 

T  * 

term  in  Z  Z on i nr 

b)  r<x/2t:  the  1/r'  terms  in  E^  and  E^,  and  the  1/r" 

term,  in  H  are  dominant. 

Since  we  know  that  power  drops  off  as  the  inverse  square  of 
distance,  and  that  the  cower  flow  at  a  field  ocint  can  be 


we  can  extrapolate  from  a)  above  (E  a/n«  k/r  ,  where  k  =  a 
constant)  that  the  far-field  begins  at  r  =  X/2ir  and  continues 
out  to  infinity,  while  the  near-field  begins  at  the  monopole 
and  continues  out  to  r  =  X  /  2  7 r.* 

It  is  interesting  to  note  the  complementary  relationship 
that  is  maintained  by  the  power  in  this  transition  region. 
Instantaneous  power  flow  at  a  point  in  space  is  given  by  the 
complex  Poynting  vector: 


3*  =  1/ 2 (E  x  H) 


For  the  short  monopole,  the  complex  Poynting  vector  is 


oriented  in  the  radial  direction  and  giver,  bv : 


rv,  _  -  j  v,  ->  Q 

( (  3  in  7)  e  '  )  “  (  — 


The  real  part  of  this  equation  gives  the  average  rate  of  flow 
of  Peal  power  (*l/r“)  while  the  Imaginary  part  is  proportions 
to  the  difference  between  average  magnetic  energv  dens it v 


ms  argument  is  0 

hi  Th 


f fared  not  as  a  proof,  but  as  a  reasonab 
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and  average  electric  energy  density  (*l/rw).  Interpretation 
of  this  equation  indicates  that  for: 

a)  r>\/2ir:  Reactive  power  flow  decreases  rapidly 

r- 

Real  power  flow  decreases  not  so  rapidly  ( < *^) 

r‘ 

b)  r  =  a/2it:  .Average  Real  power  flow  and  average  Reactive 
power  flow  are  equal 

c)  r<\/'2i:  Reactive  power  flow  becomes  much  greater 
(*r°)  than  Real  power  flow 

The  transverse  component  of  the  complex  Roynting  vector: 


cases,  the  current  distribution  in  conductors  is  complex. 

This  complexity  is  the  result  of  non-uniform  currents  and 
charge  densities,  capacitance  and  inductance  generated  by 
these  non-uniformities,  and  discontinuities  from  impressed 
forcing  functions. 

It  is  well  established  that  as  the  S.eld  observation  point 
moves  within  the  X/2-t  distance  from  the  antenna,  than  the 
measured  values  deviate  from  those  predicted  by  Equations  (5) 
through  (7).  It  is  also  well  established  that  Maxwells 
equations  are  valid  for  all  space.  Consequently,  the  con¬ 
clusion  is  that  the  sinusoidal  current  distribution  is  r.o~ 
sufficiently  accurate  for  determining  field  points  wirhin 
the  near-field  region.  The  problem  rhe.n  i;  :o  find  a  better 
description  of  the  current  distribution.  This  is  the  domain 
of  contemporary  analysis. 

Before  proceeding,  it  is  necessary  to  reflect  chat  the 
analysis  performed  thus  far  has  been  for  a  short  antenna. 

The  model  for  this  study  is  a  resonant  nor.opcle;  when  the 
antenna  is  lengthened  to  resonance  some  interesting 
characteristics  appear.  Primary  among  these  is  the  f  ac~ 
that  the  reactive  component  of  the  imoedance  vanishes ,  which 
says  the  input  reactance  goes  to  cero.  The  current  on  the 
antenna  and  the  voltage  required  to  place  it  there  are 
determined  solely  by  the  resistance,  the  real  part  of  the 
impedance.  When  a  sinusoidal  current  is  assumed  on  the 
antenna,  the  tangential  field  component  is  given  bv: 

2  1 


I 


The  (-2cosSh)  term  representing  radiation  from  the  base  of 
the  monopole  has  disaDpeared.  The  monopole  appears  as  a 
point  source  radiating  from  its  height,  and  an  image  source 
found  at  the  reflection  through  the  ground  plane. 

Also,  the  resonant  condition  occurs  at  a  length  less 
than  that  expected  from  the  X/^  relaf ionshio .  This  is  the 
result  of  an  "effective"  lengthening  of  the  antenna  by 
capacitance  and  "fringing"  at  the  end-cap  of  the  wire. 
Capacitance  is  caused  by  a  build-up  of  charge  at  the  end  of 
a  wire  of  finite  width.  The  change  in  direction  of  the 
geometry  at  the  end  of  "he  wire  also  changes  the  direction 
of  the  lines  of  force  emanating  perpendicular  from  it;  this 
is  the  fringing  affect.  These  factors  cause  resonance  to 
occur  for  an  antenna  length  of  h  =  \  -  5  where  5  is  a 

function  of  the  capacitance  caused  by  the  twc  effects 
discussed  above. 

The  short  rnonopole  development  discussed  above  is  very 
restrictive  and  cannot  be  applied  directly  to  resonant 
moncooles.  However,  when  the  resonant  peculiar  character! s 
tics  mentioned  are  taken  into  consideration,  the  resulti.or 
development  will  follow  the  same  crccedure  as  used  for  ~h'J 
short  mcncnole.  for  "he  resonant  motorola  lorment  an: 
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CONTEMPORARY  ANALYSIS 


Schelkunoff  [Ref.  3:  pp.  370-374]  identifies  three  methods 
for  obtaining  a  closer  approximation  to  the  actual  current 
distribution  on  an  antenna.  They  are: 

a)  The  Integral  Equation  method 

b)  The  Sweep-Off  method 


c)  The  Mode  Theory  of  Antennas  method. 

Each  are  discussed  below,  in  succession. 

1 .  The  Integral  Equation  Method 

The  assumption  is  made  that  the  current  is  district 
on  the  surface  of  a  hollow  cylinder  which  is  divided  into 
filaments  of  angular  density  (I(:')/2t)  and  angular  width 
( d$ ' ) .  The  Kernel  Function  is  modified  to  allow  for 
variations  in  3  and  is  expressed  as  an  integral  around  the 
circumference  of  the  wire  as: 


(’-»)  =  u_  •  ( -  -  *  ^  J,  *  )  *  f 


=  wav  a  function. 


e  field  i.ntensitv  for  a  straight  current  filament: 


(1— b  +  3  '  .)  KM)  dc’ 
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rhe  geometrical  conr iguration  or  tne  system;  hence,  a 
generalised  solution  requires  a  number  of  terms  which  ramie 
approaches  infinity.  This  overwhelming  problem  of  account! 
for  all  terms  limits  the  attractiveness  of  this  method. 

The  Integral  Equation  method  appears  most  premising 
for  generating  a  better  approximation  to  the  current.  Many 


different  aooroaches  to  solving  the  ecuati 


. eve  seer. 


prcocsed  and  virtually  all  of  them  lend  themselves  to 
solution  by  digital  commuter.  Three  of  these  aooroaches  wi 


oe  discusser 


ne  next  section. 
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'E'  field  outside  the  conductor  (expressed  in  terms  of 
vector  potential,  ’A').  The  wave  equation  is  solved  as  a  sur 
of  a  complimentary  function  and  a  particular  integral.  A 
constant,  C0,  in  the  solution  is  evaluated  in  terms  of  incut 
conditions  at  the  terminals,  and  the  vector  potential,  A,  is 
expressed  in  terms  of  the  antenna  current.  C0  and  A  are 
inserted  in  the  solution  for  the  wave  equation  obtaining  an 
Integral  Equation  in  current,  which  is  Hailed s  Equation.  A 
partial  solution  for  the  current  is  obtained  by  evaluating 
one  of  the  integrals  so  the  current  is  expressed  as  the  sum 
of  several  terms,  some  of  which  also  induce  the  current,  I. 
neglecting  certain  terms  in  I,  an  approximate  (zero  order) 
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which  he  claims  is  more  accurate,  since  continuitv  of  high 
order  derivatives  oc  '  I'  do  not  have  to  be  soecificallv 
considered  at  segment  boundaries  when  solving  the  Integral 
Equation.  He  Then  assumes  a  Parabolic  Interne lative  Fi~  fi 
the  current  distriburion  over  each  segment  (the  distribute 
over  the  whole  antenna  is  not  necessarily  parabolic)  where 

O  X"*0  C  i  3  0  V  31 2.U  0  3  o  **  C  !u  v"* -2  "tr  cl  v‘>0  ^  ^  3.C  V  ^  2  2 ,o  c- 

He  also  defines  three  moment  functions  that  describe  the 
Kernel  Function  and  its  variation  with  respect  to  i>  and  I 
over  a  segment.  The  non- zero  field  points  are  then  extras 
as  double  summations  of  the  moment  functions  over  each  int- 
vai.  It  avoid  the  double  summations,  he  s slits  the  Kernel 
Function  into  two  factors  (a  technique  used  by  Fallen,  as 
well)  which  facilitate'.-  the  commutation  of  the  dominant  te: 


/  HzMe^2’  dt’  = 


A  /  *  h  c  ( £  ) 


where:  s  (C )  is  the  Fourier  Transform  of  the  current 

n 

harmonics,  A  is  an  infinite  series  of  constants.  Ecuaticn 
n 

(2^)  is  substituted  into  Equation  (22)  and  solved  by  applying 

Talerkins  Method.*  This  yields  an  infinite  system  of 

equations  which  can  be  truncated  and  solved  for  the  unknown 

A  .  A  similar  treatment  can  be  aotlied  to  the  ether  ncn-zerc 
n  •  ■ 

field  equations  to  cast  them  into  solvable  form.  However, 
the  equations  so  defined  converge  slowly;  therefore,  addi¬ 
tional  manipulation  is  required.  The  additional  manioulation 
becomes  rather  involved  and  will  not  be  discussed  here. 


Di  scus s ion 


All  three  of  the  above  methods  ar 
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sound,  al_  three  give  reasonable  resu.ts  over  _arge  portions 
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numerical  solution  ov  oomouter. 
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amount  of  valuable  computer  resources  ( time ,  core,  and  cost) 
consumed  in  making  these  calculations.  ?aloano  indicates  the 


larkin’  3  Method  so'. 


in  Integral  Equation  via  the  Method 


’orient  s  to  be  discusses  r.el:w.  .-peciri  oa__y ,  2<a_er<in  s 
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expressions  he  derives  for  field  components  are  slow  in 
converging,  and  even  after  modifications  to  his  basic 
derivation  to  hasten  the  convergence,  he  states  that  mere 
than  1500  harmonics  (a  1500th  order  linear  system  of  equations) 
are  needed  to  satisfactorily  match  boundary  conditions. 

Chang  does  not  indicate  directly  the  amount  of 
computer  resources  consumed  by  his  method,  but  he  doe-;  imply 
his  method  is  adopted  because  a  simpler  method  (The  Method  of 


Monen t  s ) 


inaccurate : 


"...in  the  numerical  computation,  the  integral  from  -h  to  h 
is  usually  subdivided  into  a  finite  number  cf  segments;  in 
each  segment  the  current  is  interpolated  by  a  polynomial 
with  coefficients  expressible  in  terms  of  values  of  current 
at  a  few  sample  points  within  that  segment.  The  higher 
derivatives  of  I  "(a)  are  therefore  discontinuous  at  boundary 
points  between  any  two  segments."  [Ref.  5] 

Cns  of  the  intentions  of  this  study  is  to  show  that 
“he  numerical  Electromagnetics  Cede  (h'EC)  ,  which  uses  the 
Me*hod  of  Moments ,  has  been  mcdicied  to  account  rcr  the 
discontinuities  cf  the  first  derivative  at  segment  boundaries, 
and  that  this  mcdificaticn  is  sufficient  tc  viald  results  that 
are  accurate  for  --r.sine^rint  anelvsis. 


4 


structure,  an  incident  plane  wave  with  either  linear 


elliptic  polarization,  or  the  field  may  be  due  to  a  Hertzian 
dipole.  Output  may  include  current  and  charge  density,  power 
gain  or  directive  gain,  near-  or  far-zone  electric  or  magnetic 
fields,  impedance  or  admittance,  total  radiated  power  or  input 
power.  It  is  suitable  for  either  antenna  analysis  or  cross- 
section  scattering  and  electromagnetic  pulse  (Id?)  studies. 

The  code  will  handle  junctions  of  wires  of  uneven  radius,  free 
ends  of  wires  that  have  finite  radius,  wires  of  variable  radius 
and  coupling  between  wires. 

MIC  utilizes  the  dauss-Poolittie  Method  for  solving  the 
matrix  equations  generated  by  the  Method  of  Moments  when 
solving  the  integral  equations.  It  also  allows  for  use  of 
rotational  or  plane  symmetry  to  reduce  computation  time.  When 
the  impedanc-  matrix  is  too  large  to  be  contained  in  core.  MIC 
has  the  otticn  of  storing  portions  out  of  core.  It  allows  the 
' self-interaction  matrix*  for  a  structure  to  be  computed, 
factored  for  solution,  and  stored  on  tape  or  file.  A  solution 
for  a  new  antenna  that  enters  this  environment  requires  only 
one  evaluation  of  the  'self-interaction  matrix'  for  the 


n  — 
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plates)  to  which  the  Electric  Field  Integral  Equation  (EFIE) 
or  Magnetic  Field  Integral  Equation  (MFIE)  are  aoplied.  Is 
with  any  description  or  the  real  world,  there  are  aocrcxina- 
tions  ,  but  because  of  the  versatility  in  modeling  the  geo¬ 
metry  of  a  structure,  its  approximations  more  closely 
resemble  nature.  The  resemblance  is  strongly  influenced  by 
the  choice  of  coning  (i.e.  dissecting)  the  structure  in  the 
program.  The  smaller  the  geometric  elements,  the  closer  the 
model  comes  to  reality.  However,  the  smaller  the  elements, 
the  larger  the  number  of  elements,  which  means  the  larger  the 
matrix  of  equations  and  hence,  the  more  costly  the  solution. 
There  is  a  point  beyond  which  smaller  cones  do  not  yield  a 
substantially  more  accurate  solution;  it  may  be  necessary  to 
rerun  the  problem  with  increasingly  smaller  elements  to  find 
the  ocir.t  of  diminishing  returns.  The  choice  of  orooer  coning 
then  is  gained  with  experience  and  becomes  as  much  of  an  art 
as  it  is  a  science.  The  guidelines  for  the  science  aspect 
are  as  follows. 

1 .  i  re  o 

'egmerro  sncuid  follow  the  paths  of  conductors  using 
a  piece-wise  linear  fit  or.  curves.  lenerallv,  segment  lengths 
(l)  should  be  less  than  .11;  shorter  segments  (.051)  or  less 


\  depends  on  the  Kernel  used  in  the  Integral  Equation.  Two 
options  exist.  The  thin-wire  Kernel  models  a  filament  curren 
while  the  extended  Kernel  models  a  uniform  current  distribu¬ 
tion  around  the  segment  surface.  The  field  of  the  distribu¬ 
ted  current  is  approximated  by  the  first  two  terms  in  a 
series  expansion  of  the  exact  field,  in  powers  of  a".  The 
first  (a"1)  term  is  identical  to  the  thin-wire  Kernel;  the 
second  term  extends  the  accuracv  for  lamer  values  of  a. 


'action  II. A)  and  both  require  2na/\<<i.  The  thin-wire 


Kernel  requires  a  d/a>5;  the  extended  Kernel  relaxes  this 
to  d/a>2.  These  values  ensure  errors  are  less  than  1 ^  .  The 
extended  Kernel  is  used  at  free  wire  ends  and  be “ween  taralle 
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_na  naranetsr  terming  a  idr:n  is  a  norma^  unis  vector, 
originating  from  she  center  of  the  catch,  defined  in  Cartesia 
Each  catch  must  be  connected  bv  a  wire  at  the 


sates 


patches,  and  the  function  is  numerically  integrated.  fat : h 
with  wires  connected  to  then  (active  patches)  should  ce 
chosen  to  be  approximately  square  with  sides  parallel  to  th 
unit  vectors  defining  the  surface.  Inly  one  wire  may  ;c r.r.e 
to  a  patch,  that  wire  may  not  be  connected  to  another  patch 
and  it  may  not  lie  in  the  plane  of  the  patch.  A  minimum  of 
about  25  patches  should  be  used  cer  wavelength  of  surface 
area;  the  maximum  sice  for  an  individual  catch  is  about  .2- 
squara  wavelengths.  The  number  of  patches  used  increases, 
and  the  sice  of  each  patch  decreases,  as  the  radius  cf 
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restricted  to  modeling  voluminous  todies  with  closed  surfac 
parallel  surfaces  on  opposite  sides  cannot  be  too  close 
t  cgether . 

1  .  Ground  Plane 


s.  ^ -5 C  1  y u 

the  ground;  r.owev 
where  a  -  wire  ra 
ground  plane,  and 
A  finitel 

im.ase  modified  bv  one  rresne 
c sufficients .  This  method  i 


c  xose 
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having  a  large  horizontal  extent  over  the  ground.  The 
Sommerfeid/Nort on  node  1  uses  the  exact  solution  and  is 
accurate  close  to  the  ground;  the  horizontal  restriction  is 
the  sane  as  for  a  perfect  ground.  This  model  is  only  used 
for  wire-wire  interactions,  for  surfaces  the  code  reverts  to 
Fresnel  Reflection  coefficients.  At  the  present  time,  a 
ground  stake  cannot  be  used,  but  wires  nay  end  on  a  perfectly 
conducting  ground  plane  if  the  derivative  of  the  current  at 
the  end  of  the  wire  is  zero.  For  wires,  there  are  options  for 
a  radial-wire  ground  screen  anproxi.natior.  and  a  twc-nedium 
ground  approximation  based  on  modified  reflection  coefficients 
These  options  allow  considerable  savings  in  computational  tine 
when  limited  accuracy  can  be  tolerated. 

1.  INTEGRAL  EQUATIONS 

1 .  Electric  Field  Integral  location  (FFTF) 

The  FFIF  is  used  for  thin-wire  structures  ’with  snail 
conductor  volume.  I-  is  derived  from  the  electric  field 
t i *" i cn  f r. t  di ..‘inbuticn  r  o  n  ~ insd  tc 
?urf-ac9  d5  a  oerftcolv  luc  Tins'  body  : 

:(?)  =  JLj  Tjr1)  T( r ,r * )  c.V  C5' 

where:  r*  is  the  source  point,  r  is  she  ocservat^  n  point. 
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impedance  of  medium.  The  principal  value  integral  f  is 
indicated  so  r-*-r'  in  rhe  limit. 

An  integral  equation  for  the  current  induced  on  3  by 
an  incident  field  E1  can  be  obtained  from  Equation  (25)  and 
the  boundary  condition: 


i  (  r )  X  [  E  (  r )  +  E  ~  ( r  )  ]  =  3 


where  reS,  n(r)  is  the  unit  normal  vector  at  r,  E"  is  the 

field  due  to  the  induced  current  5  .  This  internal  equation 
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is  given  bv: 
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D.  METHOD  OF  MOMEMTS  AMD  THE  MATCH  POIMT 

The  Method  of  Moments  is  a  technique  whereby  an  integral 
equation  is  reduced  to  a  system  of  linear  algebraic  equations 
which  are  easily  handled  by  high  speed  digital  computers. 

1.  Mathematical  Concent 
The  method  at; 
co0r-3.tr or  z  z  ths  fo rn  i 
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function  which  is  an  unknown  response  to  an  rrxcita4-: 
■which  is  given,  and  found  in  the  range  of  1. 
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- ne  unknown  function,  r,  may  be 
M  linearly  independent  Basis  functions  sp 
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f  a  solution  to  Equation  (-3)  exists  and  is 
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,  then  r be  inverse  operator,  L  ~  ,  also 


such  that 
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(  44  ) 


which  is  a  solution  to  Equation  (43). 

The  efficiency  of  computations  and  accuracy  of 
solution  is  largely  dependent  on  the  choice  of  3asis  Function. 
Factors  which  should  guide  this  choice  are: 

a)  accuracy  of  solution  desired 

b)  ease  of  evaluation  of  matrix  elements 

c)  size  of  the  matrix  that  can  be  inverted 

i)  the  realization  of  a  'well-conditioned'  matrix. 

There  are  two  general  classes  of  Basis  Functions  from  which 
to  choose ,  entire  domain  and  sub-domain.  The  sub-domain  class 
has  fewer  elements;  hence,  its  execution  time  is  less,  i~ 
simplifies  the  evaluation  of  the  inner  product  integral  and 
still  ensures  the  matrix  Cll  will  be  'well-conditioned'. 

2  .  b'ZT's  Application 

,'iEC's  interrelation  of  current  on  wires  enforces 


(35),  or  Equation's  (36)  to  (33);  the  equation  selected 
depending  on  the  structure  being  modeled.  These  integral 
equations  are  the  basis  from  which  the  system  of  linear 
equations  will  be  generated. 

A  source  model  with  a  known  excitation  voltage  can 
be  related  tc  the  E-field  impressed  on  the  antenna;  hence, 

E~  (=  e  from  Equation  (39))  is  known.  '.ve  know  we  are  looking 
for  a  more  accurate  description  cf  the  current  distribution 
on  the  antenna;  hence,  the  current  I  (=  f  from  Equation  (39)) 
is  the  unknown  response.  The  operator  1  (an  integral 
operator  in  this  case)  is  the  composite  of  all  other  terms/ 
factors  (other  than  E  and  I)  found  in  the  integral  equations. 
Hence,  these  terms  can  be  cast  into  the  form  of  Equation  (39) 
as  : 


onto  segments  to  which  it  is  directly  connected,  going  t-: 
zero  with  zero  derivative  at  the  outer  ends  of  connected 


segments , 


Ecuation  (53)  is  then  substituted  into  Eauaticn  (43) 


to  give: 


M  M 

(z)  =  L[  I  I  g^(z')J  =  Z  I  .  L [ g . ( z ' )  ] 
i =1  3  J  j=l  J  J 


where:  I.  and  the  summation  have  been  moved  outside  *:he 

domain  of  L  as  a  result  of  its  linearity. 

The  next  task  is  to  identify  a  set  of  weighting 
functions,  {w.'.  When  the  weighting  functions  are  set  equal 
to  the  basis  functions  (i.e.  w.  =  g.),  this  is  known  as 
laierkins  Method."  In  MIC,  w.  is  chosen  as  a  set  of  Dirac 
Delta  Functions: 


W.  (z)  =  '  ( Z-  =^) 


where : 
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encroach  (lee  Iec~i:n  II. 1.1; 


on  E ■" , 

which  then  st 

ates  that  the 

tangential 

scattered  field. 

(which  generates 

the  current) , 

is  equal  to 

the  negative 

of  the 

known  impress 

ed  field,  E". 

Stated  in 

equation  form: 

-,i 

—  -.3 

(57) 

Hence  , 

the  s^/stem  of 

equations  can 

be  solved 

for  the  exact 

current 

at  the  match 

points,  and  b 

y  enforcing 

current  and 

charge 

continuity  at 

segment  ends , 

a  very  cio 

se  approximation 

to  the  actual  current  distribution  can  be  obtained. 

Equation  (54)  is  solved  by  the  Gauss -Doolittle  method. 
In  this  method  the  matrix  [E]  is  defined  as  two  sub-matrices: 

[I]  =  ['}]  (59) 

where  sub-matrix  [A]  with  elements  A.,  reoresents  the  electric 

i  - 

field  at  the  center  cf  segment  i  due  to  the  5th  basis 
function,  centered  on  segment  j.  Sub-matrix  [Sj,  with  elements 
2.. ,  reoresents  magnetic  fields  generated  by  segment  basis 
functions.  In  the  execution  of  this  method,  the  matrix  [I] 
is  factored  into  the  product  of  an  utter  triangular  matrix 
[  ■  J  3  and  a  lower  triangular  matrix  Cl],  That  is: 


l  -  J 


[Ll  C'J 3 


(50 


Equation  (£4)  then  becomes: 


[V]  =  El]  [L]  [U]  (6'3) 

and  [I]  is  computed  in  two  stages: 

[V]  =  [L]  [?1  (51) 

and : 

[F]  =  [U]  Cl]  (52) 

where  Equation  (51)  is  solved  first  for  [F]  (an  intermediate 
matrix)  bv  forward  substitution,  and  then  for  Cl]  by  back-war 
substitution. 

how  that  the  current  has  been  determined,  the  other 

non-cero  scalar  field  com.Donenfs  (E  and  H,  )  can  easilv  be 

r 


compute a 


IV.  MODEL  DESCRIPTION  AMD  COMPUTATIONAL  RESULTS 

A.  MODEL  DESCRIPTION 

The  basic  model  used  for  this  study  is  a  simple  nonoooie 
cciinear  with  the  I-axis  and  mounted  perpendicular  to  a  per¬ 
fectly  reflecting  ground  plane,  that  is  itself  located  in  the 
X-Y  plane.  As  can  be  seen  frcm  Figure  1,  the  antenna  is  75 
meters  in  height  and  .6  meters  in  diameter.  It  is  radiating 
at  a  frequency  of  I  Megahertz  (  the  middle  of  the  AM  broad¬ 
cast  band)  with  a  wavelength  (\=c/f )  of  300  meters;  it  is  a 
quarter-wavelength  moncpoie  operating  at  its  resonant  fre¬ 
quency.  Transmit  power  is  fed  to  the  model  through  its  base 
by  a  voltage  that  is  normalised  to  obtain  1003  watts  of 
radiated  tower. 

Virtually  ail  modeling  schemes  start  from  the  basic 
premise  that  I-  and  H- field  component  values  can  be  determined 
from  a  knowledge  of  the  current  distribution  on  the  antenna. 
Three  generalised  methods  for  describing  the  current  distri¬ 
bution  will  be  discussed  here:  1)  assuming  a  current 
sistr ibution ,  0)  modal  modeling,  and  •)  sub-domain  modeling. 
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approximated  by  geometric  distributions  (sin,  cos,  and  linear 
constants)  on  the  elements.  This  is  the  method  used  by  TIC 
[Ref.  2],  The  NEC  code  uses  the  scalar  form  of  the  integral 
equations  to  determine  the  current  in  each  of  the  elements, 
then  combines  the  results  from  each  substructure  to  seduce  the 
performance  of  the  total  structure. 

3.  MODEL  CONFIGURATIONS 

In  this  study  three  configurations  of  the  model  were 
assumed  and  values  for  the  tangential,  radial,  and  teak  (time 
and  spatial)  fields  were  calculated  for  each.  The  three 
configurations  were  for  a  structure  modeled  with  5,  15,  and 
25  sub-domains.  In  each  of  the  configurations,  a  null  value 
for  the  tangential  E-field  was  imposed  at  the  center  point  of 
each  sub-domain  (this  point  is  called  the  match  point).  Since 
a  closed  form  solvable  function  has  not  been  identified  which 
describes  the  current  on  the  antenna,  while  maintaining  a 
tangential  E-field  equal  to  cere,  it  has  become  necessary  to 
aooroximate  the  current  distribution.  This  is  done  in  NEC 
through  the  Method  of  Moments  as  discussed  in  Section  III.D. 
Lsing  this  method  the  current  can  be  anoroxim.atec  to  an 
arbitrary  degree  of  accuracy  by  increasing  N.  F: rturately , 
reasonable  accuracy  is  obtained  in  most  situations  with  a 
relatively  small  number  of  segments,  tvpicaliv  N  =  13/\. 
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C.  EXPECTATICNS/03SERVATI0NS./ ANOMALIES 

This  study  utilized  three  conf igurations  to  determine  two 
relevant  facts:  a)  whether  the  calculated  value  of  "he 
near-field  was  decendent  on  the  number  of  segments  used, 
and  b)  whether  the  size  of  the  segment  conformed  to  the 
guidelines  used  in  the  far-field. 

A  non-analytic  overview  of  the  slots  generated  by  this 
study  indicated  that  the  calculated  value  of  the  near-field 
was  independent  of  the  number  of  segments  into  which  the 
structure  was  divided.  This  was  determined  by  overlaying 


l 


plots  of  five  and  fifteen  segments  on  a  twenty- five  segment 
oiot.  An  example  of  one  of  these  plots  is  shown  in  Figure  u . 
The  spurious  field  emitted  from  the  base  of  the  antenna  has 
significantly  different  values.  These  variations  were  the 
result  of  using  different  values  of  input  voltage,  when  the 
configuration  was  changed,  to  maintain  a  normalized  radiated 
cower  output  of  1  300  watts."  The  values  along  -he  antenna 
and  at  its  top  are  essentially  "he  same ,  for  respective 
oositions  on  the  antenna,  for  all  three  configurations . 

The  EEC  manual  [Ter.  1]  recommends  segment  lengths  ( „) 


than  r  or  norma,  mode. mg,  an  1 


in  this  study,  segment  lengths 


330/13  =  : 3 


■am 


urn  re censed 


segment  length. 
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from  the  base  of 


nation  was  obtained  from  a  single  height  plot  that  was  not 
available  on  the  double  height  plot,  so  only  plots  to  twice 
the  antenna  height  are  discussed  in  this  section. 

Additionally,  closely  spaced  tiers  within  a  segment  wer 


2.  E-rielJ  clots  (E_,  E  .  E  ,  )  as  a  function  of  ?.  to 

c  r  peak 

a  distance  of  one  (1)  wavelength  from  the  monopole,  for 
heights  of  2/130ths,  1/4,  1/2,  1,  and  2  times  the  height 
of  the  monopole. 

3.  K-field  plots  (H)  as  a  function  of  Z  to  twice  the 
height  of  the  mcr.opole,  at  distances  of  1,  2,  £,  13,  130, 
and  1003  radii  from  the  axis. 

u.  H-field  plots  (H)  as  a  function  of  E  to  a  distance  of 
one  (1)  wavelength  for  heights  of  2/lOOths,  l/u,  1/2,  1, 
and  2  times  the  height  of  the  monetcle. 

Each  of  these  categories  are  discussed  in  further  detail 
be  low . 

1 .  E(Z)  Plots 

As  shown  in  figure  5,  £__( Z)  along  the  antenna  surface 
was  much  lower  in  magnitude  (about  60  dS  lower)  than  its  value 
at  the  top.  These  lower  values  varied  by  an  order  cf  magni- 
3.3  <3  r03w2.tr  of  !II2C  forcing  E  -  j  ^ n3*"cf'i  oc '  nt ,c: 
approximating  values  (from  the  current  distribution )  else- 
’.vP.3r0*  The  se  rdiicdl  v  3  r  i  3  tr  i o n 3  v;$p6  tuicki'?  3r 

the  observation  point  ().?.)  moved  a way  from  the  surface.  At 
a  distance  of  five  radii  (figure  6),  the  integrating  effect 
:f  a  wider  field  of  view  had  essential!-/  smoothed  out  the 


ve,  ana  auso  aecreasea  tn 


:q  H  4-’io  a  ^  ~ 


rerence  m  maenituce  oe tween 


va.jas  a* 


:he  top  of  the  antenna  and  dhcse  along  ads  length. 

?.  moves  to  a  distance  of  130  radii,  the  integrating 
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the  antenna,  and  this  uniform  magnitude  will  continue  to 

decrease  in  absolute  value  as  the  distance  of  the  0.?.  moves 

radially  from  the  antenna. 

The  magnitude  of  E_,  and  E^  are  approximately  equal 

at  the  top  of  the  monoDole,  but  E  maintains  a  value  within 

r 

an  order  of  magnitude  of  its  teak  value  all  along  the  antenna 
surface  (Figure  ").  As  a  result  of  boundary  conditions  on  E^ 
for  the  ground  olane  (similar  to  those  for  Z_  on  the  antenna 
surface),  the  magnitude  of  I  quickly  approaches  zero  at  the 
base.  The  shape  of  this  curve  remains  pretty  much  the  same 
as  observations  are  mace  at  greater  distances.  The  magnitude 
of  the  curve  decreases  slowly  and  approaches  a  uniform  value 
along  the  antenna. 

E  is  a  cornoosife  of  E  and  E,.  Tear  the  antenna, 

pea.<  -  r  z 

clots  of  E  ,  (Figure  3)  strong!*/  resemble  clots  of  E  .  Thi 
teak  -  •  r 

implies  that  E  is  the  dominant  ccm.Dcner.t  in.  the  near-riels 

re-*  ion  •  However,  as  the  3.?.  '  s  are  moved  further  cutwarc, 

the  E  .  /E  resemblance  -ades  and  there  is  a  ccint  around 
o  <2  cL<  r 

\ /Em  (see  Figure  3)  where  the  dominance  begins  to  shift  to 
Z_ .  Then  E„  becomes  completely  dominant,  the  value  ot  E^ea- 
is  accroaching  a  uniform  value  along  the  antenna.  Ehss 
results  from  the  electromagnetic  waves,  generated  by  the 
ooint  source  at  the  top  of  the  antenna,  approaching  a  plane 
wave  front  at  large  distances.  Also  noted  in  Figure  3  is  the 
sharp  decline  of  E  above  one  height  or  the  antenna.  -his 
ir  dilates  very  little  radiation  in  the  axial  direction,  t.oe 
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greatest  percentage  is  directed  radially  outward.  A  nore 
complete  set  of  Ed)  Plots  are  found  in  Figure  A-l  through 
A-lS  in  Appendix  A. 

2 .  Z(R)  Plots 

Figure  10  shows  the  very  large  value  of  E^  found  at 
the  top  of  the  antenna  decreases  rapidly  in  magnitude  as  the 
0. P.  moves  radially  away  from  the  monopole .  At  a  radial 
distance  of  approximately  X/2-rthe  rapid  decline  in  Z_^  sub¬ 
sides  and  approximates  the  attenuation  with  distance 
observed  at  other  points  along  the  antenna  (see  Figure  11). 
Since  E„  =  C  at  the  surface,  as  shown  in  Figure  11,  the 
magnitude  of  the  field  rises  rapidly  with  radial  distance  to 
a  peak  at  about  20  meters  from  the  antenna,  from  which  the 
magnitude  slowly  declines. 

The  rate  of  decline  of  E^(R)  appears  to  be  uniform 
for  all  positions  along  the  antenna,  the  only  difference 
being  the  magnitude  from  which  the  decline  begins.  (Figure  12 

Close  to  the  antenna,  the  E  .  curve  follows  the  E 

pea.<  r 

curve.  At  seme  ^r.kTiOwn  point,  E„(R)  has  diminished  in 

magnitude  less  than  E  ( R) ;  hence,  the  E  .  curve  switches  to 
-  r  oea.t 

follow  E„(R)  (Figure  1  ’  > .  The  point  at  which  this  switch 
occurs  falls  at  a  greater  radial  distance  for  2. P.'s  at 
greater  heights.  This  is  the  result  of  a  larger  starting 
magnitude  (at  treater  heights)  from  which  E  basins  its 
j r. i.  ^  OT'm  v'1 -ti *r  7* or'  1  ^  ~ 
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Plot  Er  v»  r  (or  a  modoi  of  15  logmontt. 

H  of  1/2  th#  mono  polo  kotgkt  along  tko  r  axtt. 


larger  radial  distance  is  allowed  for  it  to  decline  before 

the  magnitude  of  E,  exceeds  E  .  The  point  at  which  this 

occurs  is  so  sharoly  defined,  that  the  curve  of  E  .  has  an 

■  J  ’  peaK 

"apparent  discontinuity"  in  it. 

At  twice  the  monopole  height,  the  curve  of  Eq 
(Figure  14)  has  an  inverted  oust.  This  results  from  the  same 
shift  in  E  /E  dominance  discussed  above.  Directly  above  the 


antenna,  E  =  0 :  hence,  E  is  larger  than  -  .  But  as  the  3. 

’  '  '  'T  -0 


moves  radially,  E_  ieclir.es  v/hil- 


rises  to  a  oea.c. 


forms  the  inverted  cuso:  then,  E  also  begins  a  decline,  but 
net  as  raoid  as  E  .  More  clots  of  E(E)  are  found  in  Figures 
A-19  through  A- 3 3  in  Appendix  A. 

3.  H(Z)  Plots 

At  the  surface  of  the  monopole,  the  magnitude  of  Hv 

r- 

is  virtually  constant;  there  is  a  slight  decrease  near  the 
top  (Figure  15).  Above  the  monooole,  there  is  a  share  drop 


(about  30  i5 ’ s )  in  intensity 


y  or  r.  ;  t.ne 


field  is  yirtuall” 


r.on-exis  tent 


:imil  a: 


"eater  sis~ances 


from  the  antenna  follow  the  same  tvoe  of  curye ,  extent,  the 
magnitude  of  the  fielo  oerallel  to  the  antenna  decreases  with 


mcreas  mm 


ranees  from  the  antenna,  and  the  sharp  decline 


in  intensity  above  the  height  of  the  antenna  becomes  less 
abrupt.  At  a  distance  of  one  wavelength,  the  difference 


oam.e. 


and  ■ibov® 
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H ( R)  Plots 


4  . 

Again,  at  the  surface  of  the  monopole  the  magnitude 
of  H  is  virtually  constant.  The  magnitude  of  H  decreases 
rapidly  with  increasing  radial  distance,  and  uniformly  for 
various  antenna  heights.  Figure  16  is  an  example  curve  which 
delineates  the  behavior  of  H  all  along  the  axis.  As  indicated 
above,  for  O.P.'s  directly  above  the  antenna,  H  is  virtually 
non-existent.  As  the  O.P.  moves  radially  outward,  at  twice 
the  height  of  the  antenna,  K  rises  quickly  to  a  peak,  then 
slowly  declines  in  magnitude.  More  plots  of  H ( R )  are 
found  in  Figure  A-u3  through  A-^4  in  Appendix  A. 


A  VALIDATION  OF  NEC  NEAR-FIELD  CALCULATIONS 


NEC  is  still  evolving.  On-going  efforts  toward  refining 
and  improving  its  calculations  are  constantly  occurring.  In 
its  present  form,  it  has  been  extensively  used,  and  docu¬ 
mented  to  provide  accurate  results  in  the  far-field;  its 
validity  in  close  to  the  antenna  is  not  so  well  documented. 
To  show  that  NEC  provides  reasonable  results  in  near-field 
regions,  three  tests  have  been  run.  They  are  described  in 
the  following  paragraphs. 

A.  COMPARISON  OF  NEC  WITH  CLASSICAL  ANALYSIS  RESULTS 

Mills  [Ref.  3],  using  equations  similar  to  those  found 
in  Ref.  4  ( pp .  323  to  333),  calculated  theoretical  values 
for  as  a  function  of  radial  distance  from  the  antenna, 
for  snort  antennas  of  various  lengths.  Two  of  his  plots  are 
of  interest  to  this  study;  they  are  for  diooles  of  half- 
length  (h)  equal  to  .13  and  .33.  The  results  of  his 
calculations  are  clotted  as  broken  lines  in  Figure  17,  the 
clot  of  NEC’s  calculations  for  the  same  antennas  are  shown 
as  solid  lines.  Remembering,  the  assumption  of  sinusoidal 
current  is  accurate  onlv  for  the  quarter-wavelength  mcnopole 


:VAs  the  antenna  becomes  shorter,  the  sinusoidal  peak  occurs 
further  back  on  the  feed  line  and  the  distribution  on  the 
antenna  becomes  more  linear. 


it  is  interesting  to  note  NEC’s  calculations  are  more  closely 

matched  as  the  antenna  length  aoproaches  the  resonant 

condition,  especially  for  radial  distances  greater  than 

X/2tt  (=  300/2  =  48  meters  for  this  study).  Balzano  [Ref.  7] 

points  out,  that  inside  this  distance,  the  classical  approach 

breaks  down,  "...because  nowhere  in  the  near  field  (p<X/2-rr) 

does  the  dipole  look  like  an  elementary  source  (infinitesimal 

dimension)  from  the  observation  point".  Balzano  also  points 

out  that  the  classical  approach  has  a  constantly  increasing 

E„  field  ud  to  the  dioole  surface  which  is  inconeruent  with 
z. 

energy  density  boundary  conditions  at  the  antenna  edge. 
Therefore,  the  leveling  off  of  the  magnitude  of  E„,  as 
oredicted  by  NEC,  is  the  more  reasonable  plot. 

3.  COMP ARISEN  0?  NEC  WITH  EMPIRICAL  ANALYSIS  RESULTS 

Balzano  in  Ref.  7  (Part  II)  performed  some  measurements 
on  a  dipole  constructed  from  bronze  rods  of  .3015  meters 

_3 

diameter  (=  3.75  X  13  \)  that  -was  .133  meters  in  length 

(=  .i712\),  including  a  .3018  meter  gap  at  the  center  from 
which  it  was  fed. 

A  model  of  this  antenna  was  constructed  using  NEC.  The 
configuration  chosen  was  a  moncpoie  of  19  segments,  each 
segment  being  .30^9  meters  long.  The  monopole  was  mounted 
abov"  a  perfectly  conducting  ground  plane  in  a  manner 
similar  to  that  used  for  this  study.  Plots  of  E_  and  E^  were 
run  as  a  rune  tier,  cf  E  for  radial  distances  of  .  33  3  ,  .  3  36  , 
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.01,  .02,  .03,  .04,  .05,  and  .06  meters  from  the  dipole  axis. 

The  results  of  these  plots  are  shown  in  Appendix  B.  With 
the  exception  of  a  scaling  factor,  the  information  contained 
in  these  plots  is  the  same  as  that  discussed  in  Section  IV. C. 
for  the  study.  To  further  compare  the  results  of  NEC  with 
the  measured  values,  the  measured  peaks  of  the  square  of  the 
electric  field  intensity,  for  both  axial  and  radial  polariza¬ 
tions,  at  the  radial  distances  indicated  above,  were  extracted 
from  plots  in  Ref.  7.  These  values  were  plotted  in  Figures 
13  and  19  along  with  the  square  of  the  peak  of  the  fields 
extracted  from  the  plots  of  the  model  run  by  NEC.  As  can  be 
seen,  very  close  correlation  between  the  model  and  the 
measurements  were  obtained.  The  largest  discrepancy  was  at 
.003  meters  distance  where  the  calculated  value  was  inordinately 
higher  than  the  measured  value,  relative  to  other  calculated/ 
measured  pairs.  Balzano’s  calculated  value  was  also  quite  a 
bit  higher,  which  would  indicate  the  measured  value  was  a  bit 
low.  Credence  for  this  argument  is  supported  by  Balzano's 
statement,  "...at  closer  distances,  capacitive  effects 
became  predominant". 

:.  COMPARISON  Or  NEC  WITH  MAXWELL' S  EQUATIONS 

Chang  [Ref.  5 ]  ncints  out  that  "...on  the  surface  of  the 
mor.onole  the  field  component  H ,  should  be  proportional  to  the 
current  distribution  I(Z )...".  The  NIC  calculated  value  of 
H ,  was  ronoared  to  the  current  flowing  on  the  model.  The 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


From  the  results  of  Section  V,  it  can  be  seen  that  NEC  is 
an  effective  analysis  tool  for  near-field  calculations  ana 
can  be  used  and  trusted  to  give  meaningful  results.  The 
tool  aannot  be  applied  thoughtlessly ;  care  must  be  exercised 
to  ensure  the  model  accurately  represents  th^  actual  physical 
structure.  Zoning  considerations  must  be  factored  into  the 
procedure  to  ensure  that  sufficiently  small  segments  are  used 
to  get  the  desired  resolution,  but  not  so  many  as  to  be 
wasteful  of  resources.  Electrical  parameters  must  be  con¬ 
sidered  to  ensure  they  represent  the  operational  environment 
into  which  the  structure  will  be  immersed.  Then  these  items 
are  adequately  considered,  NEC  will  provide  the  desired 
results  in  a  timely  and  cost  effective  manner. 

Areas  of  further  study  could  include  expanded  investi¬ 
gations  into  the  maximum  segment  length  for  near-field 
calculations ,  investigation  of  the  characteristics  of  the 
field  emitted  from  the  base  insulator  region,  and  effects 
on  the  near-field  of  using  multiple  parallel  wires  with  cross 
connections  in  the  manner  in  which  broadcast  towers  are 
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